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Abstract 

Using 24 million ip 1 = ip(2S) decays in CLEO-c, we have searched for higher multipole admixtures 
in electric-dipole-dominated radiative transitions in charmonia. We find good agreement between 
our data and theoretical predictions for magnetic quadrupole (M2) amplitudes in the transitions 
ip' ~^ 7Xci,c2 and Xci,c2 — ► jJ/4>, m striking contrast to some previous measurements. Let b% 
and a~2 denote the normalized M2 amplitudes in the respective aforementioned decays, where the 
superscript J refers to the angular momentum of the XcJ- By performing unbinned maximum 
likelihood fits to full five-parameter angular distributions, we found the following values of M2 
admixtures for J x = l: a J 2 = l = (-6.26 ± 0.63 ± 0.24) x 10~ 2 and b$ =1 = (2.76 ± 0.73 ± 0.23) x 10~ 2 , 
which agree well with theoretical expectations for a vanishing anomalous magnetic moment of the 
charm quark. For J x = 2, if we fix the electric octupole (E3) amplitudes to zero as theory predicts 
for transitions between charmonium S states and P states, we find a{= 2 = (-9.3 ± 1.6 ± 0.3) x 10~ 2 
and bi= 2 = (1.0 ± 1.3 ± 0.3) x 10~ 2 . If we allow for E3 amplitudes we find, with a four-parameter 
fit, a{= 2 = (-7.9±1.9±0.3) x 10~ 2 , b{= 2 = (0.2±1.4±0.4) x 10~ 2 , a^= 2 = (1.7±1.4±0.3) x 10~ 2 , 
and bi= 2 = (-0.8 ± 1.2 ± 0.2) x 10~ 2 . We determine the ratios a^ =1 /a^= 2 = 0.67±g;i| and 
a 2 =1 /^2 =1 = ~~ 2.27^0 99) where the theoretical predictions are independent of the charmed quark 
magnetic moment and are ai =l /ai= 2 = 0.676 ± 0.071 and a^ =1 /b( =1 = -2.27 ± 0.16. 
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I. INTRODUCTION 



The radiative transitions between spin-triplet charmonium states are known to be dom- 
inated by electric dipole (El) amplitudes, but higher multipole contributions, magnetic 
quadrupole and electric octupole (M2 and E3), are sometimes allowed. These higher multi- 
poles give information about the magnetic moment of the charm quark. To search for these 
contributions, we studied the radiative decay sequences 

e+e" ^7* -> i/j' = ip(2S) 

i>' ~> 7 ; X(cl,c2) 
X(cl,c2) -» jJ/lp 

J/ij) — > e + e~ or 

using the helicity formalism developed in Refs. [l|, 0, H, HJ . As shown in Fig. [TJ the particles 
i 3 ' i X(ci,c2), an d J/ij) are the 2 3 Si, l 3 P(i i2 ), and 1 3 Si charmonium states, respectively. For 
the J x = 1 decay sequence, we search for two multipole amplitudes, b{ =1 and a^ =1 , which 
are respectively the M2 amplitudes for the ip' - > l'Xci (b for before the Xc) an d Xd ~ > iJ/i* 
decay (a for after the Xc)- Similarly, for the J x = 2 decay sequence, we search for two M2 
amplitudes (b^ 2 , ai =2 ) an -d two E3 amplitudes (b^ =2 , ai =2 )- 
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FIG. 1: Charmonium energy levels. Only the transitions studied in this article are shown. 

The multipole amplitudes are calculated from a maximum likelihood fit of the joint 
angular distribution of the two photons 7' and 7, described by five angles for each event. 
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Polar and azimuthal angles (9', 0') denote the direction of the initial e + e~ axis relative to 
the first photon 7' (in the if)' frame), an angle # 7 y describes the direction between the two 
photons (in the Xc frame), and polar and azimuthal angles (9,(f>) denote the direction of 
the final lepton pair (£ + £~) axis relative to the second photon 7 (in the J /if) frame). These 
angles are illustrated in Fig. [2J 
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FIG. 2: Reference frames defining the angles used in this analysis. In the if)' frame, the angles 9', (f)' 
are the polar and azimuthal angles of the beam pipe (specifically, the positron's direction) relative 
to 7' defining the z'-axis, and 7 lying in the x'-z' plane (with a positive x'-component). In the 
Xc frame, the angle # 7 y is the angle between the two photons. In the J /if) frame, the angles 9, (f> 
are the polar and azimuthal angles of the two leptons (specifically, the positive lepton's direction) 
relative to 7 defining the z-axis, and 7' lying in the x-z plane (with a negative x-component). 

In previous experimental studies of Xcj — > jJ/ifi, the magnetic quadrupole amplitude in 
the decay sequences involving Xci was found to be consistent with zero, while that found 
via Xc2 was found to be several standard deviations from zero. However, theory predicts the 
ratio of these two magnetic quadrupole amplitudes to be of order unity. With CLEO's large 
sample of if)' decays, the question is ripe for re-investigation. The present paper describes 
that effort. 

Section [TT] sets the theoretical stage for the investigation. Prior experimental results are 
reviewed in Sec. IIHI The CLEO detector, data sets, and Monte Carlo samples are described 
in Sec. [TV] Sec. [V] discusses selection criteria, while Sec. |VT] is devoted to fits to the data. 
Systematic uncertainties are treated in Sec. IVII[ while Sec. IVIIII concludes. 



II. THEORETICAL CONTEXT 
A. Allowed radiative transitions 

For the radiative decays between a 3 Si state and a 3 Pi state, only El and M2 transitions 
are allowed. For 3 Si — > 3 P2 transitions, from conservation of angular momentum and parity, 
we would expect that the £"3 transition would be allowed, but this is forbidden under the 
single-quark radiation hypothesis @, 0|. Single-quark radiative transitions must have |AS| < 
1 and parity-changing transitions must have |AL| = 1, so the photon cannot carry off three 
units of angular momentum [lj]. However, for the J x = 2 case, electric octupole transitions 
are allowed if either S state has a small D admixture p], or if the P state has a small F 
admixture. There is evidence 0, 0, 0] that the if)' state is actually a mixture cosy? |2 3 Si) — 
sin ip \ 3 Di) with <p = (12 ± 2)°, so we may expect a small 63 transition amplitude. 
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B. Joint angular distribution 



The formalism developed in Refs. [l|, 0, 0] is used to construct the joint angular distribu- 
tion of the decay sequence. We denote the signal decay as 

V/(A') - 7 V) + XV) (1) 
X (u) -> 7 (a0 + WA) (2) 

with helicities in parentheses and the helicities associated with the if)' decay labeled with 
primes. For the if)' (xc) decay sequence, the helicity amplitudes are labeled B v i (A v ) and 
the multipole amplitudes are labeled bj i (aj ,), where J 7 is the angular momentum carried 
by the photon 7. The helicity amplitudes are specified by only one helicity, since parity 
conservation allows the independent helicity amplitudes to be defined for J x > v > as 

B v > = B v ^i = {—l) Jx B- u i_i , 
A v = A v , x = {-1) J *A- Vi - X . 

Here the second index of the two-index helicities refers to the photon. To form the joint 
angular distribution the if)' and J / if> density matrices must be constructed from the directions 
of the two electrons forming the if)' and the two leptons that decay from the J /if). 1 

The angles 9', <$' contain information on the polarization of the if)', while 9, <f> contain 
information on the polarization of the J /if). The angle 7 y, defined by the angle between the 
two photons in the \ c rest frame, gives information on the necessary rotation between the 
two reference frames. Frames for construction of these five angles have been shown above 
in Fig. [2J The joint angular distribution is therefore 

W (cos 9' , 4> 1 cos 7 y , cos 9, <f>) oc 

is'v'; n' = ±1 
vv\\x = ±1 



where d J J u are standard Wigner (^-functions 10 



In e + e~ — > 7* — > ip' , the polarization of the ip' along the beam axis is ±1, so the density matrix giving 
the polarizations in the direction of the beam axis (the z-axis) is given by p( x x ^ = e^ X ^e[ X ' + "'e^ \ 
where is the polarization vector for helicity A defined with components — (— 1, — i, 0)/V2, 
e (°) = (0,0,1) and e^ 1 ) = — eW* = (1,— i,0)/s/2. Generalizing to an arbitrary direction n = 
(sin 9' cos 4>', sin 9' sin (f>', cos 0'), we find that the density matrix p for ip' is 

p^'' x '\e', 4>') = £ < (V) ef > ($« - nV) . 
Similarly for the J/ij) with rh = (sin 9 cos <p, sin 9 sin cos 9), we find the density matrix is 

pi^(e, 4>) = £ 4 (X) 4 X) ( sij ~ mW ) • 
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The helicity amplitudes A v (with < v < J x ) are related to the multipole amplitudes 
aj 1 (with 1 < J 7 < J x + 1), using the Clebsch-Gordan coefficients (ji,mi, j 2 ,m2\J, M) , by 



(4) 



This expression leads to the following relationships for the J x — 1 and J x = 2 cases, 
respectively: 



AJ=l 

A( =1 



( At 2 



A J=2 

A{= 2 




(5) 



(6) 



The relationships between B v i and bj , are identical; just swap all A v and aj i with B v i and 
bj f/ in Eqs. fll])-©. These transformation matrices are norm-preserving, since the matrices 
are orthogonal. 

C. Quark magnetic moments 



If we define El, M2, and E3 to be the electric dipole, magnetic quadrupole, and electric 
octupole amplitudes, respectively, the magnetic quadrupole amplitudes are related to the 
anomalous magnetic moment of the charm quark k c by 



,j=i 



,J=2 



b J=i 



b J = 2 



M2 



En, 



VE1 2 + M2 2 4m c 
M2 



[1 + K c 



3 E~ 



VE1 2 + M2 2 + E3 2 
M2 



V / 5 4m c 



'l + K c 



<l + K c 



y/El 2 + M2 2 4m 
M2 

VE1 2 + M2 2 + E3 2 ~ V54:m c 



3 Ey 



1 + K c 



(7) 
(8) 
(9) 
(10) 



These expressions are correct to first order in i? 7 /m c or Ey /m c , assuming that the ip(lS,2S) 
are pure S states (no mi xing with D states) and that the Xc states are pure P states (no 
mixing with F states) 0, Til]. 2 

These first order relationships are derived from the non-relativistic interaction Hamilto- 
nian for photon emission from a +2/3 charged quark: 



2m, 



■(A*-p + p- A*)-fjur-K* 



(11) 



Note the misprint in [1 11 ] for their equation (41) describing a'{ to first order. This misprint was previously 
noted in footnote 1 of Ref. 12 1. 
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where e c = ||e|, \i = (e c /2m c )(l + k c ), A* is the vector potential of the emitted photon, and 
H* = V x A* is the magnetic field of the emitted photon (both A* and H* are complex 
conjugated since the photon is outgoing). We omit a spin-orbit term of the same order which 
does not affect the M2 contribution [2], 13| . 

The ratios of the predicted multipole amplitudes given by Eqs. 
of m c and k c to first order: 



at 1 ' 




E J , =1 V5 
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0.676 ±0.071 



-2.27 ±0.16 



1.000 ±0.015 



-0.297 ±0.025 



- (TTOjl are independent 

(12) 
(13) 
(14) 
(15) 



As the individual amplitudes have corrections of order (E^/m c ) 2 , we conservatively assigned 
each multipole amplitude a fractional uncertainty equal to (E 1 /m c ) 2 (using m c = 1.5 GeV, 
k c = 0) which was the dominant source of uncertainty in Eqs. (|T2|) - (|T5|) . 

The E3 amplitudes are expected to be small in view of the few-percent admixture of the 



l 3 Di state in the ip'. Although they are found to be complex in Ref. [ll|, we shall include 
them in fits assuming that they are real. 



D. Lattice QCD predictions 



Dudek et al. [141 . [15J performed lattice QCD calculations for the charmonium radiative 
transitions X(ci,c2) "~ ¥ lJ/4>- They ran lattice simulations at various values of Q 2 (the square 
of the four- vector of the photon, which is for real photons) and extrapolated to Q 2 — > 0. 

For the transition Xci — ► J J/ ip, when extrapolating the El and M2 amplitudes to Q 2 — » 
individually, they found that 

M2(Q 2 -^ 0) -0.020 ±0.017 . . 

El^Q 2 -> 0) ~ 0.23 ±0.03 —° ° 9± ° ° 7 ("J 

They concluded that data points at smaller Q 2 and improved knowledge of form factors were 
needed to make a meaningful comparison with experimental values [15j. Similarly for Xc2 
jJ/if), they found the normalized multipole amplitudes behaving as at 2 ~^ —0.39 ± 0.07, 
af 2 -> 0.010 ± 0.011 as Q 2 -> 0. 



III. PRIOR EXPERIMENTAL RESULTS 

Tables HI and ITU summarize the results from previous experiments for J x = 1 and J x = 2, 
respectively. 

For transitions involving Xci, the Crystal Ball experiment at SPEAR used 921 events of 
e + e~ — > tj)' — > j'xd —> l'lJ/4> ~^ to measure both multipole amplitudes. The 
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E-835 experiment used 2090 pp — > Xci ~~ * iJ/i* ~~ * 7e + e~ events to measure the multipole 
amplitude a J 2 =x . The CLEO-c data sample has ~40000 e + e~ — > ^i' — > ■y'xci — > j'jJ/ip 
r y' r y£ + £~ events after applying selection criteria. 

For transitions involving Xc2, the Crystal Ball experiment used 441 events of e + e~ — > ip' — > 
•y'Xc2 — > r y' r yJ/ip — >■ r yy'£ + £~ to measure both multipole amplitudes. The E-760 and E-835 
experiments used 1904 and 5908 — > x C 2 — > ^fJ/tp — > 7e + e~ events, respectively, to measure 
the multipole amplitude a^ 2 - The BESII experiment searched for the multipole amplitudes 
in a novel method looking at 418 ij}' —>■ 7Xc2 —> 77r + 7r~ events and 303 ip' — > 7Xc2 - ► ^K + K~ 
events, and the BESII fit also found a value of bi= 2 = -0.0271^9 • The CLEO-c data sample 
has ~20000 e + e~ —> ip' — ► 7'xc2 — ► I'lJ/^ — > 7 / 7C~£~ events after applying selection 
criteria. 

TABLE I: Previous experimental values vs. theoretical predictions for the normalized magnetic 
quadrupole amplitude for the decays Xci — ► lJ ( a 2 =1 ) an d i>' — ► l'Xci (^2 =1 )- 

Experiment «2 =1 H =1 Signal Events 

Crystal Ball [16, 17] -0.002^;^ 0.077lg;^ 921 

E-835 [12] 0.002 ± 0.032 ± 0.004 2090 

Theory (m c = 1.5 GeV) -0.065(1 + k c ) 0.029(1 + k c ) 



TABLE II: Previous experimental values vs. theoretical predictions for the normalized magnetic 
quadrupole amplitudes for the decays Xc2 — ► lJ/4 1 i a i =2 ) an d ifi' — ► l'Xc2 iH =2 )- 

Experiment «2 =2 ^2 =2 Signal Events 

Crystal Ball [16, 17] -0.3331^ 0.132lg;g?f 441 

E-760 [18] -0.14 ±0.06 1904 

E-835 [12] — 0.093±g;g|? ± 0.006 5908 

BESII [19] -0.05lig:oie 731 

Theory (m c = 1.5 GeV) -0.096(1 + k c ) 0.029(1 + k c ) 



Many of these experimental results disagreed with theory which predicted ratios given in 
Eqs. ffT2l- f|T5l) . The ratios of the averages of previous experimental values compared with 
theory values are 
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The first two ratios which involve the multipole amplitudes that have the most statistical 
significance strongly disagree with their theoretical predictions. As the ratios are indepen- 
dent of m c , k c and any specific quarkonium potential model to first order in £' 7 /(4m c ), we 
expect good agreement between theory and experiment. 



IV. DETECTOR, DATA, AND MONTE CARLO 
A. The CLEO detector 

Data were acquired at the ip' resonance at ^Js = 3.686 GeV using the CLEO-c detector 



located at the Cornell Electron Storage Ring (CESR), a symmetrical e + e~ collider [2d. 121 
The solid angle for detecting both charged and neutral particles is 93% of 4ir. The pho- 
tons were detected as showers in a Csl (Tl) calorimeter consisting of 7784 crystals, which 
achieved a photon energy resolution of 2.2% at 1 GeV and 5% at 100 MeV. The azimuthal 
and polar angular resolution for 100 MeV photons is cr^, azim ~ 11 mrad (19 mrad) and 
cr0 polar ~ 0.8a^ azim sin ^poia,. (10 mrad) in the barrel (endcap) region of the crystal calorime- 
ter. Charged particles were detected using a set of two cylindrical drift chambers enclosed 
within a superconducting solenoid with a 1.0 T magnetic field directed along the beam axis. 
The outer drift chamber achieved a momentum resolution of ~0.6% at p = 1 GeV and an 
azimuthal and polar angular resolution of 0"0 azim m 1 mrad and erg ^ ~ 4 mrad 20| . (In 
this article, c = 1 in mass and momentum units.) The inner six-layer stereo drift chamber 
is used to accurately measure the location of charged particles along the beam axis. 



B. Data sets and expected number of events 

For our analysis, we used the recent CLEO-c data set taken at the ip' events consisting 
of a sample of (24.45 ± 0.49) x 10 6 ip' events with a total luminosity of 48.07/pb (22|. Using 
known branching fractions T(| and the known sizes of the CLEO data sample, we can expect 
that 91900 ± 6600 J x — 1 signal events and 48200 ± 3600 J x = 2 signal events are originally 
present in the data sample. 



C. Phase space Monte Carlo 

For each of the decay sequences (J x = 1 and J x = 2), a 4.5 million event phase space 
Monte Carlo (MC) data sample was generated. The phase space MC was generated with 
EvtGen [23[ with final state radiation simulated with PHOTOS (24T |. 

The purposes of the phase space Monte Carlo are threefold. First, it is used to account 
for the variable angular efficiency of the detector after selection criteria have been applied, 
when performing the maximum likelihood fit (see Sec. I VI Al) . Second, the phase space MC 
events are used to simulate signal MC with non-zero multipole amplitudes, 02, hi (and 03, 63 
for J x = 2) via the rejection method. This is achieved by taking the five angles 9', <p', 77 /, 8, <f) 
for each phase space event and calculating the probability of that event occurring at those 
angles for the PDF W (fi; Aq) with the input values of the multipole amplitudes Aq. The 
probability for the event occurring at that angle is then compared to a random number 
uniformly distributed between and 1. Then, our simulated signal MC obeying the PDF 



9 



W(Q; Aq) consists of the events that are more probable than the corresponding random 
numbers. The third purpose of the phase space Monte Carlo is to generate projections to 
overlay upon histograms of data values. For example, after a fit to data extracts values of 
a 2 ,b 2 for a J x = 1 fit, the phase space MC can be used to generate projections in the five 
angles with the fitted values of a 2 , b 2 to be compared with the data. 



D. Generic Monte Carlo 



In order to properly simulate feed-across into the selected data sample from non-signal 
ip' decays, a "generic" MC sample was prepared. This sample consists of approximately 
120 million ?y decays, using our best estimate for all measured branching fractions 0,0, 



25l . l26l . 1271 . l28l . l29l . l30l . l3ll . |32j, |33j and matrix elements for the deca ys of ijj' and its decay 
products; unmeasured hadronic decays are simulated using JETSET 34] . The signal events 
(■?// — * 7'x(ci,c2) — ► I'lJ/^) were replaced with phase space MC events selected to have the 
desired a 2 and b 2 admixture (via the rejection method as described in Sec. IIVCI) . 



V. SELECTION CRITERIA 

Tuning of the selection criteria was designed to eliminate non-signal "impure" background 
events, while selecting the largest number of signal events. For kinematic regions in which it 
was uncertain how to apply selection criteria, we attempted to minimize the quadrature sum 
of the statistical uncertainty from signal events and the systematic uncertainty from impure 
events. Many of the starting points for our selection criteria are taken from a CLEO-c study 
22| of ip' — > h + J/ip branching fractions that included our signal decays. 



All tracks and showers investigated are required to pass standard CLEO-c criteria prior 
to any attempts at kinematic fitting. For tracks, we ensure that the track originated from 
near the interaction point (tq < 2 cm and \zq — Zi. p .| < 10 cm), is from the well-modeled 
region of the barrel (| cos 6^0^! < 0.83) or the endcap (0.85 < | cos 6^0^! < 0.93), and has 
a momentum between 1% (18.4MeV/c) and 120% (2.21GeV/c) of the beam momentum. 
The requirement for a shower is that it is not matched to a track, has | cos 6^0^-! < 0.79 or 
0.85 < I cos^poiarl < 0.93, and has an energy between 1% and 120% of the beam energy. 

All candidate events require at least two tracks and two showers to be identified. The 
two tracks and two showers used (if more are present) will be those with the greatest en- 
ergies. Two kinematic fits are then performed to generate the four 4-vectors used in the 
analysis. First, a 1C kinematic fit to the J/ip mass is performed starting with the two 
tracks, allowing shower(s) identified as bremsstrahlung photons to be associated with a 
track. Bremsstrahlung photons are identified if a shower that is not matched to a track is 
located within 100 mrad of the initial momentum vector of a track. If bremsstrahlung pho- 
tons are identified, the lepton four- vector used is the sum of the kinematically fit 4-vectors of 
the lepton plus all associated bremsstrahlung photons. Second, a 4C kinematic fit to the ip' 
4- vector is performed and the result of this fit is then subjected once more to the original 1C 
fit. The ip' 4- vector is calculated from the angle at which the electron and positron beams 
intersect (4 mrad) and the beam energy of the given run. For both the 1C and 4C kinematic 
fits, we require the reduced y 2 for both the vertex and kinematic fit to be less than 16 as 
shown in Fig. [31 This value was found by minimizing the quadrature sum of the impurity 
systematic uncertainty and statistical uncertainty. 
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To identify signal events through the J x radiative cascade, we require the reconstructed 
Xcj mass to be within 15 MeV of the true Xcj mass as constructed by adding the J/ip and 
7 four-vectors together: 

m xcj = \/\Pj/tP+P~i\ 2 = \J\Pe+ + Pe~ +P 7 I 2 ( 21 ) 

We do not apply a selection criterion based on the other Xcj mass reconstructed from the 
ip' and 7', as the 4C kinematic fit ensures that this criterion is redundant (see Fig. H]). 

Signal events must also have the J/ip decay to e + e" or ^ \T , so we require the two tracks 
to be well-identified as both being electrons or muons. We achieve this by looking at the ratio 
of the energy deposited in the calorimeter to the momentum of the track (E/p). We identify 
both tracks as electrons if the smaller E/p ratio is greater than 0.5 and the larger E/p ratio 
is greater than 0.85. Similarly, we identify both tracks as muons if (^/p) sm aiier < 0.25 and 
(E/p)\ avger < 0.5. This results in a clean e-/i separation. 

To restrict major sources of background, we apply additional criteria to address the modes 
with large branching fractions: 

B{ip' -»• ttV J/ip) = (16.84 ±0.33)% 
B(ip' — »• r]J/ip) = (3.16 ±0.07)% 
B(ip' -> 7T°J/ip) = (1.26 ± 0.13) x 10~ 3 (22) 

The dominant background mode ip' — > 7T 7T J /if) — > 7777^ + £ is reduced by requir- 
ing the third most energetic shower in the event (excluding those photons identified as 
bremsstrahlung photons) to have an energy of less than 30 MeV. To reduce the contributions 
of the background modes with monochromatic J/ip momentum, ip' — > rjJ/ip — > 77 J/ip and 
ip' -> n° J/ip -> 77J/^ where p(J / ip%, = 199 MeV and p{J/ip%,^j^ = 528 MeV, 
we require the J/ip momentum to lie between 240 MeV and 510 MeV. Note that the signal 
transition generates no events with a J/ip momentum below 238 MeV (318 MeV for J x = 2) 
or above 542 MeV. 



VI. FITTING THE DATA 

A. Basic approach and procedure 

We find the multipole amplitudes by performing a maximum likelihood fit of the selected 
data events to the probability distribution function (PDF) W(Q; A) given by Eq. (jSJ). Events 
are selected according to the criteria described in Sec. [V] and each event is described by a set 
of five angles Q = (6', <p', # 77 , 9, (p) defined in Fig. H The PDF W{Q; A) gives the probability 
for an event with angles Q to occur given a set of multipole amplitudes A = (a{,bj). The 
PDF in Eq. ([3]) is written in terms of helicity amplitudes, but can be written in terms of 
multipole amplitudes as W(Q;A) using Eq. (j3J). The total likelihood for Nd data events to 
be described by W(Q; A) is 

N d 

C w (A) = Y\W(n d ;A). (23) 

d=l 

The initially unknown angular detector efficiency e(fl) describes the probability that an 
event occurring at the angles fl will be registered by the detector and pass the selection 
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FIG. 3: Maximum reduced \ 2 in & U kinematic fits (including vertex fits) in generic Monte Carlo and 
data. Events with a maximum reduced x 2 below 16 (the dashed vertical line) are kept. Cumulative 
totals for the number of signal and impurity background events are also plotted for each potential 
value of a maximum reduced \ 2 - ( a ) Jx = 1 an< ^ O 3 ) Jx = 2- 
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FIG. 4: Plot of the Xcj mass as calculated from subtracting the four- vector of the 7' from the ijj' 
four- vector. This variable is not used as a selection criterion because the 1C and 4C kinematic fits 
ensure that this criterion is redundant with the Xci mass selection criterion generated by adding 
the J/ij) and 7 four- vectors, (a) J x = 1 and (b) J x = 2. 
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criteria. We define a new normalized PDF to account for this detector efficiency e(Q): 

F(a-A)= (24) 

and note that the original PDF VT(fi; A) is of the form 

W(fl; A) = £ cuajhkGijuin). (25) 

ijkl 

The functions G^kiiP) are obtained from the expression for W(Q;A), so this form allows 
separation of the parameters being determined (the multipole amplitudes A) and the data 
points (the angles Q). This allows us to write the denominator of the PDF in Eq. (1241) as 

i(Q')W(Q';A)dn' = J e(n / )53a i a i 6 fc 6,Gfy JH (n , )tZn' 

ijkl 

= Y^ a i a jhk J e{Q!)G ijk i{Q!)dQ! 

ijkl 

= ^ aidjbkbilijki 

ijkl 

where the detector-efficiency-dependent integrals I^i = J e(Q')Gijki(Q')dQ' are independent 
of the fitting parameters A. The integrals Iijki can be approximated by a Monte Carlo 
numerical integration technique. Using a large sample of phase space Monte Carlo events 
(Sec. IIV CI) generated uniformly in the five angles (cos#', </>', cos0 7 y, cos 9, 4>), we record 
whether each phase space MC event is reconstructed and passes the selection criteria. Using 
the known angular functions Gij^Q), we approximate the integral as 



1 T ijkl 



= J e(n')G m (Q')dn' 



N, 



phsp 



J2 ®{p)G ijkl {%) (26) 



^Vphsp p= i 

where Q(p) is 1 (0) if the pth phase space event is (not) reconstructed and iVph sp is the total 
number of phase space events. 

The most likely form of the parameters A given the PDF F(Q; A) is found by maximizing 
the logarithm of the likelihood, which is given by Eq. f[2Uj) with the PDF F instead of W. 
The logarithm of the likelihood that the parameters A in the PDF F(Q; A) describe the 
data events occurring at angles Qd is 

N d N d 

\og£(A) = io g n^;^) = E 1 °g ir (^^) 

d=l d=l 

N d r 

= £ log e(tt d ) + log W(Q d ; A) -\og^a i( i jbkkhjM . (27) 

d=l ijkl 

The first term in log £ is independent of the A, so the log likelihood only depends on the 
detector efficiency through the phase space integrals. We reduce the number of parameters in 
the fit by recognizing that the multipole amplitudes are normalized (e.g., a\ + a\ + a\ — 1). 
This method of performing an unbinned maximum likelihood over an angularly varying 
detector efficiency was first developed in Ref. 35[. The multi-dimensional optimization of 
log £'(A) was achieved using the Minuit MlGRAD variable-metric fitting routine [36| . 
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B. Statistical results of five-angle fits 



1. J x = 1 fits 

The result of the two-parameter fit to the J x = 1 data is a 2 =1 = —0.0611 ± 0.0063, 
b 2 =1 — 0.0281 ± 0.0073, based on 39363 events. The efficiency integrals in the denominator 
were calculated by simulating 4.5 million phase space MC events taking account of the 
detector geometry and selection criteria; 39.6% of events were reconstructed. Contours are 
shown in Fig. Ufa) of \/2A log£, where A log £ is the difference in log likelihood between 
the fitted values of a2, b 2 and any other values. For a pure El transition (02 = 62 = 0) the 
value of xei = v^ATogX is 11.1. 

The projections of the data in each of the five angles may be compared with curves 
based on a pure El distribution and on the fitted M2/E1 admixture. The angle 9 is of 
particular importance as it is the angle that most clearly shows the preference of the data 
for an M2/E1 admixture over a pure El transition. The projection for cos 8' also shows 
slightly better agreement with data with the fitted M2/E1 admixture than with a pure El 
transition. For the 50-bin histograms in cos#, the reduced x 2 (x /iVd.o.f.) 3 comparing the 
data with the projection at the fitted values is 42.7/47 = 0.91, while data and the pure El 
projection have a xV^d.o.f. of 108.5/49 = 2.21. 

Using the parity transformations described in Ref . [TtJ , we can fold four of the five angles 
into the positive domain without modifying the value of the likelihood calculated through 
VF(f2; A). In Fig. EJ^a) we show that the data are well-matched with the projection in | cos^l 
with the fitted values of A, but poorly matched with the pure- .El | cos#| projection. 

When we fix the ratio of the parameters to the theoretical ratio, given by Eq. (1131) . 
°2 =1 /^2 =1 = —2.274, we can perform a one-parameter fit to the five-angle J x = 1 data set. 
The result of this one-parameter fit is a J 2 =1 = -0.0615 ± 0.0055, b J 2 =1 = -a J 2 =1 /2.27A = 
0.0271 ± 0.0024, with a value of xei = V2A log£ = 11.1) nearly identical to the results of 
the two-parameter fit. The results of these two fits are compared in Table [TTT1 

TABLE III: J x = 1 five-angle fit results. The fits were performed on 39363 events satisfying the 
selection criteria described in Sec. |V} xei = y/2A log C is the number of standard deviations by 
which the fitted value differs from the pure El value. 



Fit 


a( =1 


bt 1 


Xei 




(10- 2 ) 


(10- 2 ) 




Two- p ar amet er 


-6.11 ±0.63 


2.81 ±0.73 


11.1 


One-param. (02/62 = —2.274) 


-6.15 ±0.55 


2.71 ±0.24 


11.1 


Theory (m c = 1.5 GeV) 


-6.5(1 + k c ) 


2.9(1 + k c ) 





3 The number of degrees of freedom iVd.o.f . = ^Vbins — ^Vparams — 1 where iVbins is the number of bins in the 
histogram, and -/V params is the number of free parameters in the fit. The minus one accounts for the fact 
that the projections are normalized to contain the same number of events as the original data set. 
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FIG. 5: (a) J x = 1 and (b) J x = 2 log likelihood contours as functions of (02,62) for two- 
parameter fits. The fitted values (the solid squares) are (02,62) = (—0.0611,0.0281) for J x = 1 and 
(02, 62, 0,3, 63) = (—0.093,0.010,0,0) for J x = 2. These are, respectively, 11. la and 6.2a from pure 
El (the solid circles). The theoretical values to first order in E^/m c with k c = are given by the 
dashed lines. 
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FIG. 6: (a) J x = 1 and (b) J x = 2 projections of cos# after using parity transformations to fold 
the data set into positive cos 0', 4>' , cos # 7 y , cos 9. For J x = 1 (J x = 2) the values of X 2 /A^j. .f. 
for the 25 bin histogram describing the data to correspond with the two-parameter (02,62) fitted 
projection are 16.2/22 = 0.74 (20.3/22 = 0.92) and to correspond with the pure El projection 
are 80.29/24 = 3.35 (35.5/24 = 1.48). The fitted and pure El projections are selected from the 
same phase space MC sample (via the rejection method), resulting in the correlation of statistical 
fluctuations in the two projections. 
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2. J x = 2 fits 

As the J x = 2 PDF is parameterized by four multipole amplitudes (02, 62, «3, ^3), there 
are several choices of fits to be performed. The simplest would be a two-parameter fit with 
a 3 — ^3 — 0, as the E3 amplitudes should be zero in the absence of significant S-D state 
mixing. For this type of fit to the 19755 signal events, we find a 2 =2 = —0.093 ± 0.016, 
b J 2 = 2 = 0.010±0.013 with these fit values favored by 6.2a over a fit with pure El. 

Allowing for S-D mixing in the ip' state, the b% =2 amplitude may be non-zero. When we 
perform a three-parameter fit (setting a^ =2 = 0), we find a 2 =2 = —0.093 ± 0.016, b 2 =2 = 
0.007 ± 0.013, b(= 2 = -0.008 ± 0.011, favored by 6.3a over pure El. 

When we allow a non-zero b% =2 amplitude, but fix the ratio of a 2 =2 jb 2 =2 = —3.367 by 
Eq. (1151) . we can perform a two-parameter fit that allows for S — D mixing in the ip' state. 
The results of this two-parameter fit are a J 2 = 2 = -0.092 ± 0.016, b J 2 = 2 = -a 2 /=2 /3.367 = 
0.027 ± 0.005, b(= 2 = -0.001 ± 0.011, favored by 6.1a over pure El. 

When we perform the fit for the full four parameters (a 2 , b 2 , a 3 , 63), we find a 2 =2 = 
-0.079 ± 0.019, aj= 2 = 0.002 ± 0.014, b{= 2 = 0.017 ± 0.014, b^= 2 = -0.008 ± 0.012, favored 
by 6.4a over pure El. 

For the five-angle fit with two parameters, we plot the data with the pure El projection 
and the fitted value projection of | cos^l in Fig. [6](b). As for J x = 1, the fitted values match 
the data better than the pure El projection. 

The results of the above fits are summarized in Table HVl In all cases there is at least 6.1a 
evidence for multipoles beyond El in the transition ip' — > 7'x C 2 — y I'lJ/ip- The contours 
for ^2A log £ for 02 vs b 2 for the two-parameter fits are shown in Fig. [5](b); the contours for 
all other pairs of variables for all the fits are Gaussian-shaped with a single local maximum. 

TABLE IV: J x = 2 five-angle fit results. The fits were performed on the 19755 signal events 
satisfying the selection criteria described in Sec. |V] xei = \/2A log C is the number of standard 
deviations by which the fitted value differs from the pure El value. 



Fit 


4= 2 


bt 2 


4= 2 


H= 2 


Xei 




(10- 2 ) 


(10- 2 ) 


(icr 2 ) 


(10- 2 ) 





Two-parameter -9.3 ±1.6 1.0 ± 1.3 6.2 

Three-parameter -9.3 ±1.6 0.7 ±1.3 -0.8 ± 1.2 6.3 

Two-param. (b 2 = f§fe) -9.2 ± 1.6 2.7 ±0.5 -0.1 ± 1.1 6.1 

Four-parameter -7.9 ± 1.9 0.2 ± 1.4 1.7 ± 1.4 -0.8 ± 1.2 6.4 

Theory (m c = 1.5 GeV) -9.6(1 + k c ) 2.9(1 ± k c ) Model dep. 



VII. SYSTEMATIC UNCERTAINTIES 

We now present the results of systematic studies for the fits to the five-angle distributions 
performed in the previous section. For J x = 1, we perform all systematic studies on the two- 
parameter fit (a 2 , b 2 ), as the one-parameter fixed-ratio fit produces nearly identical results. 
However, for J x = 2, there are four types of five-angle fits: 



18 



• Two-parameter fit (02, 62) with 03 = 63 = (no S-D or P-F mixing), 

• Three-parameter fit (02, ^3) with CI3 = 0, relevant for a D-wave admixture in the ip f , 

• Fixed-ratio two-parameter fit (a 2 , 63) with 6 2 = — a 2 /3.367 and a 3 = 0, and 

• Four-parameter fit (02,62,03,03)- 



In this paper, we describe in detail the systematic studies for the J x = 1 case and the J x = 2 
case where 03 = 63 = 0. Systematic studies for the other three J x = 2 cases are discussed in 



detail in Ref. 37 



For many investigations into a possible systematic uncertainty, we perform an ensemble of 
fits on samples of signal events selected from a phase space data set via the rejection method 
to follow W(Q; Ao) for a given input set of multipole parameters. For each multipole a, we 
calculate the following parameters from the results of these ensembles of fits, with N cns MC 
events in each member of the ensemble: 

• (a), the mean of the fitted multipole amplitude over the ensemble of tests, with a 
statistical error corresponding to the variation of the fitted multipole amplitude over 
the ensemble of tests, 

• erf*, the (mean of the) nominal uncertainty from each individual likelihood fit to mul- 
tipole amplitude, 

• A/ a \, the deviation of the mean from the MC-generated value of the amplitude in units 
of the expected deviation of the mean a^ a \ = a &t / \/N ens defined as: 



(a) - a Input 



A (a) = J,, 7Tr > alld ( 28 ) 



A ff y, the deviation of the standard deviation when a potential systematic effect is 
present compared to the standard deviation without the effect being present (in units of 
the expected fluctuation the best estimate of the standard deviation from an ensemble 
of iV measurements o a = o j\JlN ( [lol ] Sec. 32.1.1), defined as: 

^j-with syst 0. without syst 

A ^ = „/ /97v— • (29) 



For a we list the mean of the nominal uncertainty, but for all tests performed the nominal 
uncertainty from every likelihood fit in the ensemble was essentially constant to the level of 
precision quoted. 

For all of the systematic tests from ensembles of measurements, we assign a systematic 
uncertainty if either (a) we find that there is a significant bias |A/ \| > 1 or if (b) there is an 
uncertainty that widens the ensemble distribution above the expected statistical fluctuation 
evidenced by A cr ( a ) > 1. 
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A. Toy MC check of fitting procedure 



To test the accuracy of the fitting procedure described in Sec. IVI Al an ensemble of toy 
Monte Carlo fitting trials was performed. For each trial, we generated a large number of 
phase space events, where each event is described by five random numbers for each of the 
variables (cos 8', <f)', cos# 7 y, cos 8, 0) uniformly distributed over their ranges. 

We generated a set of toy signal Monte Carlo events by selecting events from a separate 
MC phase space data set via the rejection method, so the events are described by W(ft; Aq) 
for an input set of multipole parameters Aq. 

To test the J x = 1 ( J x = 2) fits, we performed an ensemble of 200 toy MC trials in which 
each trial had iV S ig = 40000 (20000) signal events after selection criteria were applied. We 
analytically calculated the phase space integrals, as the toy MC was thrown at 100% detector 
efficiency. In assigning the systematic uncertainty, we set the multipole amplitudes of the 
toy signal Monte Carlo to be a( =l = -0.065, h J = x = 0.029, (a(= 2 = -0.096, b(= 2 = 0.029, 
a 3 =2 = 0; H =2 = 0) ^he expected value if k c = to first order in E 7 /m c . Fits with 
other values of input parameters recover the input results to similar precision. We find 
no systematic bias or uncertainty is associated with the fitting procedure described in this 
method, as the ensemble of trials is Gaussian-distributed with a width according to the 
statistical uncertainty. 

B. Amount of phase space Monte Carlo needed for efficiency integrals 

Using too few phase space Monte Carlo events would give poor approximations to the effi- 
ciency integrals, introducing an overall systematic uncertainty to the results of the maximum 
likelihood fit. We use 4.5 million phase space events for the normalization, approximately 
100 times the J x — 1 data set (the larger of the two). By varying the size of the MC sample, 
we find no systematic uncertainty associated with any number of events exceeding 10 5 , and 
hence assign no systematic error to this source. 

C. Impurity systematic uncertainties 

For the J x — 1 (J x = 2) selection criteria, approximately 0.23% (0.29%) of the events 
that pass the selection criteria are not signal events, but a background mode that must be 
considered for the possibility of introducing a systematic bias or uncertainty to our result. 
Taking our five-fold generic Monte Carlo data set and splitting it into five independent data 
sets, we find a purity and efficiency of 99.77% and 39.6% (99.71% and 36.0%). The main 
sources of impurity background modes for J x — 1 are if}' — > 7r°7r° J/ 'if} and if>' — > "fxd (where 
the Xci decay was not to jJ/if> followed by J /if) — > For J x = 2 they are ip' — > j'xci 

and if) 1 -> 7T°7T°J/ip. 

For each of the five independent generic MC impurity backgrounds, we perform 31 (37) 
trials with and without the impurity background events present. For each trial, we replace 
the signal events originally present with phase space events selected via the rejection method 
to come up with many independent data sets. For each trial we perform one fit with no 
impurities present and one fit with the impurities. For a given set of impure background 
events, we find that the bias due to impurities varies very little between different trials. In 
Table IVl (Table IVl)) . we list the difference from the fit with no impurities. For J x = 1, we find 
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a significant impurity bias that is relatively constant among all five sets of impure events, 
so we correct our fitted result for this impurity bias and assign a systematic uncertainty of 
half of the bias. For the J x = 2 case, we find that the impurity bias significantly fluctuates 
between background data sets, so we assign a systematic uncertainty of the size of the 
fluctuation of the impurity bias. 

TABLE V: Generic MC tests for a systematic bias from impure events for J x = 1. We split the 
five-fold generic MC data set into five data sets labeled (A)-(E), replacing the generic Monte Carlo 
signal events with events selected to obey W(fi; Aq) from the 4.5M event phase space MC data 
set. For each of these five data sets, we performed an ensemble of thirty-one fits. The difference 
rows show the shift in values of a 2 and b 2 comparing the individual fits before and after impurities 
are added. A positive shift means that to obtain the pure results we should subtract the bias from 
impurities. The set (A-E) is the result from adding all five data samples of impure events to a 
regular-sized set of signal events, and demonstrates how the impurities scale linearly in the J x = 1 
case. 



Type 


<a 2 > 

(10- 2 ) 


(10- 2 ) 


A(a 2 ) (62) 

(10- 2 ) 


(10- 2 ) 




Pure 


-6.54 ±0.50 


0.63 


-0.32 2.97 ±0.71 


0.73 


0.52 


Difference w/impurities (A) 


0.150 ±0.002 




0.058 ± 0.003 






Difference w/impurities (B) 


0.120 ±0.002 




0.053 ± 0.003 






Difference w/impurities (C) 


0.140 ±0.003 




0.060 ± 0.005 






Difference w/impurities (D) 


0.216 ± 0.004 




0.095 ± 0.005 






Difference w/impurities (E) 


0.109 ±0.002 




-0.031 ± 0.003 






Difference w/impurities (A-E) 


0.730 ±0.011 




0.241 ±0.019 






Input 


-6.50 




2.90 






(Impurity bias) 


0.15 ±0.03 




0.05 ± 0.03 







D. Final state radiation 

Another possible source of systematic uncertainty is the effect of final state radiation 
(FSR), which can alter the directions of the two leptons in the J ftp rest frame affecting the 
variables cos 6 and <fi. Generation of Monte Carlo samples has been done using EvtGen, 
which models final state radiation in the decay sequences J/ip — > £ + £~ with PHOTOS. We 
estimate the effect of final state radiation by performing signal fits on the angles Q from 
generator level four-vectors, both before and after final state radiation has been added. We 
use the rejection method to select events, so that the pre-FSR generator level four-vectors 
follow the PDF W(Q; Ao) for an input value of the multipole amplitudes Aq. We also 
use the pre-FSR four-vectors when selecting the phase space events to be used as 'signal' 
described by the PDF W(Q;A ) with a given A = (a 2 ,b 2 ) = (-0.065,0.029) (for J x = 2, 
Aq = (a 2 ,b 2 ,a s ,b s ) = (—0.096,0.029,0.0,0.0)). We then compare the fit on the selected 
events using the pre-FSR and post-FSR generator level to check for a systematic uncertainty 
from final state radiation. Comparing A CT ( a ) and A( ) for each multipole parameter for FSR, 



21 



TABLE VI: Generic MC tests for a systematic bias from impure events for J x = 2 for two-parameter 
(02,62) fit with 03 = 63 = 0. We find that the impurities add a negligible systematic uncertainty 
when compared with the statistical uncertainty. 



Tvnp 




(do) 


\ u a,2/ 


A/ \ (bo) 


(oP) A/a \ 
\ u b 2 l ^<&2> 








(10" 2 ) 


(lO^ 2 ) 


(10~ 2 ) 


Pure 




-9.8 ± 1.4 


1.6 


-0.6 3.0 ±1.3 


1.2 0.6 


Difference w/impurities 


(A) 


-0.005 ± 0.006 




0.078 ± 0.003 




Difference w/impurities 


(B) 


0.080 ± 0.004 




-0.011 ±0.005 




Difference w/impurities 


(C) 


-0.008 ±0.011 




0.149 ±0.004 




Difference w/impurities 


(D) 


0.022 ± 0.003 




-0.050 ± 0.003 




Difference w/impurities 


(E) 


-0.041 ±0.002 




0.027 ±0.003 




Difference w/impurities 


(A-E) 


0.047 ±0.019 




0.190 ±0.011 




Input 




-9.6 




2.9 




(Impurity Bias) 




0.009 ± 0.040 




0.038 ± 0.070 





we find no statistically significant evidence for a systematic uncertainty due to Final State 
Radiation. 

The angular distribution was also fit using data from J/tj) — > e + e~ and J/ij) — > fi + fi~ only, 
selected by the E/p selection criterion. Without correcting for possible systematic biases, 
we found the results in Table I VII I Preliminary studies of simulated signal data (generated 
from phase space MC) indicated that the most accurate result is obtained by performing a 
fit to the combined J/ip — ► and J/ip — > e + e~ dataset, while still blinded to the actual 
dataset. The results from the fit to the muon-only dataset were similar to the results from 
the combined dataset. The electron-only dataset produces similar results to the fit results of 
the combined dataset, except for the Xci case where the two-parameter electron-only result 
deviates from the combined result by approximately 1.4a. However, fixing the 02/^2 ratio 
reduces the deviation of the electron-only result to less than la even in this worst case. 
Therefore, we assign no additional systematic uncertainty due to FSR from the results of 
the muon-only and electron-only fits. 

E. Choice of kinematic fits 

For our final analysis, we perform a 1C kinematic fit to the J/ip mass and a 4C kinematic 
fit to the ip' four momentum of the lab frame, and also perform bremsstrahlung reconstruc- 
tion on each track if any showers were tagged as bremsstrahlung radiation belonging to the 
track. To test for possible systematic effects, we perform an ensemble of tests on phase 
space MC shaped to have A = (-0.065, 0.029) for J x = 1 and A = (-0.096, 0.029, 0, 0) for 
J x = 2 with four- vectors selected to have the pre-FSR generator photons follow W(Q; Ao). 
We construct the four- vectors for the variables in three ways: (1) Post-FSR generator level 
four- vectors; (2) 1C and 4C kinematic fits without bremsstrahlung recovery; (3) 1C and 4C 
kinematic fits with bremsstrahlung recovery. For each four-vector type, we perform as many 
fits as possible using a data size (after selection criteria) of 40000 J x = 1 (20000 J x = 2) 
events in each fit. We find no statistically significant systematic effect from this procedure. 
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TABLE VII: Fits to the angular distribution using only events where the J/ip decays to two muons 
or two electrons. The Xci (Xc2) dataset contained 20968 (10563) muon-only events and 18395 
(9192) electron-only events. 



Fit 

± lb 




rt C1 

U2 


hn 
L>2 


ft n 

"3 


ho 






\ 1U ) 


/in- 2 ! 

(IV ) 


no _2 "i 
i ±u ) 


no _2 "i 
i ±u ) 


Xcl 


± WU-pdl dill UiLl 


7 oq _i_ n 87 

— / .iO X u.ul 


o o _i_ i n 

L.L I l.U 






Xcl 


Fixed-ratio fifi 


— D.oO ± U. (0 


O.U It U.O 






Xcl 


Two-param ee 


-4.85 ± 0.92 


3.5 ± 1.0 






Xcl 


Fixed-ratio ee 


-5.36 ±0.81 


2.4 ±0.4 






Xc2 


Two-param fj.fi 


-8.1 ±2.1 


1.2 ± 1.7 








Xc2 


Three-param fifi 


-8.1 ±2.1 


1.1 ± 1.9 





-0.3 ± 1.6 


Xc2 


Fixed-ratio fifi 


-8.1 ±2.1 


2.4 ±0.6 





0.2 ± 1.5 


Xc2 


Four-param fifi 


-5.4 ±3.0 


0.0 ±2.0 


3.6 ± 1.8 


-0.2 ± 1.7 


Xc2 


Two-param ee 


-10.7 ±2.3 


0.8 ± 1.8 








Xc2 


Three-param ee 


-10.7 ±2.3 


0.2 ±2.0 





-1.4 ± 1.8 


Xc2 


Fixed-ratio ee 


-10.5 ±2.3 


3.1 ±0.7 





-0.4 ± 1.6 


Xc2 


Four-param ee 


-11.2 ±3.0 


0.4 ±2.1 


-0.6 ±2.2 


-1.4 ± 1.8 



F. Variation of selection criteria 

To look for an additional systematic uncertainty from possible variations of selection 
criteria, we looked at effects of the following variations on statistical and systematic impurity 
uncertainties: maximum third shower energy, maximum reduced \ 2 i Xc mass window, and 
maximum cosine of polar angle for photons in the barrel region. Variations were explored 
which loosened and tightened all our selection criteria. For J x = 1 we found that the default 
criteria (defined in Sec. [Vj had the smallest quadrature sum of the statistical uncertainty 
with impurity systematic uncertainty. We further found that over the ensemble of tests 
involving various criteria, the mean from the ensemble of tests for 02 and 62 (when no 
impurities were present) varied only slightly. For the J x = 2 two-parameter (02,62) fit 
case, we found that while we were quite near the minimal total quadrature sum for the 
default criteria, we could have achieved a ~ 3% improvement if we loosened these conditions. 
However, to achieve that ~ 3% improvement requires increasing the number of impure events 
by a factor of approximately five as shown in Fig. [21 so this was not performed. 

After looking at the effect of variations of selection criteria on an ensemble of tests using 
the "signal" data selected from phase space MC via the rejection method, we looked at the 
actual effect of performing fits to data after applying various criteria. These results show the 
sensitivity of the data to the chosen criteria. For the J x — 1 case shown in Table IVIIIl we 
perform the fits using the various criteria, and then correct for the impurity bias. We then 
consider the ensemble of bias-corrected data fits aed assign a systematic uncertainty using 
the standard deviation of the fitted results over the 7 types of criteria considered. We find 
a systematic uncertainty of (0.19,0.22) x 10~ 2 for (o2 =1 ,&2 =1 ) m performing fits to data. 

For J x = 2 (Table HX|) . we follow a similar procedure but do not correct for impurity 
biases before calculating the systematic uncertainty, as the impurity bias in all cases is less 
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TABLE VIII: Results of data fits for J x = 1 when applying various selection criteria. For all 
selection criteria considered, a systematic uncertainty is found of (0.19,0.22) x 10 -2 for (02,62); 
respectively, over the variation of the criteria. 



Criteria 



Default 

£3rdShwr < 18MeV 
£3rdShwr < 50MeV 

Xtf. < 10 

Xlf. < 30 

Xc mass ±10 MeV 

Xc mass ±20 MeV 

|cosCb r , r ptl<0.77 

| cos <£$|< 0.80 

Ensemble 



„bias cor 

(io- 2 ) 

-6.26 ±0.63 ±0.15 
-6.43 ± 0.64 ± 0.08 
-5.73 ±0.60 ±0.30 
-6.23 ±0.65 ±0.05 
-6.30 ±0.61 ±0.32 
-6.36 ±0.65 ±0.11 
-6.10 ±0.62 ±0.18 
-6.18 ±0.65 ±0.16 
-6.17 ±0.62 ±0.16 
-6.20 ±0.19 



^bias cor 

(io- 2 ) 

2.76 ±0.73 ±0.06 
2.67 ± 0.73 ± 0.06 
2.45 ±0.72 ±0.13 
2.33 ± 0.75 ± 0.03 
3.10 ±0.71 ±0.04 
2.85 ± 0.75 ± 0.04 
2.78 ± 0.69 ± 0.09 
2.97 ±0.75 ±0.07 
2.73 ±0.72 ±0.07 
2.74 ±0.22 



than 1/10 the statistical uncertainty, so any correction would be of very little significance. 
We find in this case systematic uncertainties of (0.3, 0.3) x 10~ 2 for (02, b 2 ) when performing 
the two-parameter fit with (a 3 = 6 3 = 0). 

TABLE IX: Results of data fits when applying various selection criteria to J x = 2 two-parameter 
(02, 62) fits (03 = 63 = 0). For all sets of criteria, a systematic uncertainty is found of (0.3, 0.3) x 10~ 2 
for (02,62)1 respectively, over the variation of the criteria. 



Criteria 


02 


62 




(IO" 2 ) 


(IO" 2 ) 


Default 


-9.3 ± 1.6 


1.0 ±1.3 


£3rdShwr < 18MeV 


-9.4 ± 1.6 


0.6 ± 1.3 


£3rdShwr < 50MeV 


-9.8 ±1.6 


0.5 ± 1.3 


xlt < 10 


-9.1 ± 1.6 


1.3 ± 1.3 


Xlt < 30 


-9.5 ± 1.5 


0.4 ± 1.2 


Xc mass ±10 MeV 


-8.7 ± 1.6 


1.0 ± 1.3 


Xc mass ±20 MeV 


-9.8 ± 1.5 


0.8 ±1.3 


1 cos Cb^K 0.77 


-9.6 ± 1.6 


1.2 ±1.3 


| cos Cb r ,prJ< 0.80 


-9.5 ± 1.5 


1.3 ± 1.3 


Ensemble 


-9.4 ±0.3 


0.9 ±0.3 



G. Summary of systematic uncertainties and biases 

The systematic uncertainties and biases for J x = 1 are summarized in Table |2D We find 
the total systematic uncertainty to be (0.24,0.23) x 10~ 2 for (a^ 1 , b^ =1 ), respectively. The 



24 



systematic uncertainties for the J x = 2 two-parameter fit (02,62) axe summarized in Table 
Kit and for other fits in Tables IXHHXIVl In each case the total systematic error is the 
quadrature sum of systematic uncertainties, and the statistical uncertainty for the data fits 
is given for comparison. We do not find any systematic biases for the J x = 2 case. 

TABLE X: Systematic uncertainties and biases for J x = 1. The total systematic error is the 
quadrature sum of systematic uncertainties and the signed sum of systematic biases. The statistical 
uncertainty from the data fits is given for comparison. 



Systematic uncertainty 


at 1 




bt 1 






Uncertainty 


Bias 


Uncertainty 


Bias 




(10~ 2 ) ' 


(10- 2 ) 


(10~ 2 ) ' 


(10- 2 ) 


Generic MC impurities 


0.15 


0.15 


0.05 


0.05 


Selection criteria 


0.19 




0.22 




Total systematic uncert. 


0.24 


0.15 


0.23 


0.05 


Statistical uncertainty 


0.63 




0.73 





TABLE XI: Systematic uncertainties for J x = 2 two-parameter fit with 02,62- 



Systematic uncertainty 


4= 2 


bt 2 




(10- 2 ) 


(10- 2 ) 


Generic MC impurities 


0.04 


0.07 


Selection criteria 


0.33 


0.33 


Total systematic uncert. 


0.3 


0.3 


Statistical uncertainty 


1.6 


1.3 



TABLE XII: Systematic uncertainties for J x = 2 three-parameter fit for 02,62,63. 



Systematic uncertainty 


a J 2 = 2 


b J 2 = 2 


bi= 2 




(10- 2 ) 


(10- 2 ) 


(10- 2 ) 


Generic MC impurities 


0.04 


0.07 


0.03 


Selection criteria 


0.33 


0.34 


0.20 


Total systematic uncert. 


0.3 


0.3 


0.2 


Statistical uncertainty 


1.6 


1.4 


1.2 
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TABLE XIII: Systematic uncertainties for J x = 2 two-parameter fit for 02,63 with fixed values of 
62 = —a2/3. 367 and 03 = 0. 



Systematic uncertainty 


a J 2 = 2 


H= 2 




(10- 2 ) 


(1(T 2 ) 


Generic MC impurities 


0.04 


0.04 


Selection criteria 


0.34 


0.23 


Total systematic uncert. 


0.3 


0.2 


Statistical uncertainty 


1.6 


1.1 



TABLE XIV: Systematic uncertainties for J x = 2 four-parameter fit with 02, 62, 03, 63. 



Systematic uncertainty 


a J 2 = 2 


bt 2 


a{= 2 


bt 2 




(10- 2 ) 


(10- 2 ) 


(10- 2 ) 


(10- 2 ) 


Generic MC impurities 


0.06 


0.08 


0.08 


0.03 


Selection criteria 


0.24 


0.39 


0.28 


0.20 


Total systematic uncert. 


0.3 


0.4 


0.3 


0.2 


Statistical uncertainty 


1.9 


1.5 


1.4 


1.2 



VIII. CONCLUSIONS 

A. Normalized multipole amplitudes 

The results of our bias-corrected fits with systematic uncertainties for J x — 1 with the 
two-parameter fit are: 

a J = x = (-6.26 ± 0.63 ± 0.24) x 10^ 2 , (30) 
b{ =1 = (2.76 ±0.73 ±0.23) x 10~ 2 . (31) 

The results of our fits with systematic uncertainties for J x = 2 with the two-parameter 
fit (a 2 , 62) with (Z3 = 63 = are: 

a J 2 = 2 = (-9.3 ± 1.6 ± 0.3) x 10~ 2 , (32) 
b J 2 = 2 = (1.0 ± 1.3 ± 0.3) x 10~ 2 ; (33) 

for the three-parameter fit (a 2 , b 2 , 63) with a 3 = 0: 

a 3 = 2 = (-9.3 ± 1.6 ± 0.3) x 10" 2 , (34) 
b J = 2 = (0.7 ± 1.4 ± 0.3) x 10" 2 , (35) 
bl= 2 = (-0.8 ± 1.2 ± 0.2) x 10" 2 ; (36) 
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for the two-parameter fit (a<z, 63) with fixed values of b 2 = —02/3. 367 and 03 = 0: 

a J = 2 = (-9.2 ± 1.6 ± 0.3) x 10" 2 , (37) 

b J=* = = (2.7 ±0.5 ±0.1) x 10~ 2 , (38) 

2 3 367 V ) V ) 

b J = 2 = (-0.1 ± 1.1 ± 0.2) x 10" 2 ; (39) 

and for the four-parameter fit (ci2, 62, «3, ^3): 

a J 2 = 2 = (-7.9 ± 1.9 ± 0.3) x 10" 2 , (40) 

b{ =2 = (0.2 ± 1.5 ± 0.4) x 10~ 2 , (41) 

al= 2 = (1.7 ± 1.4 ±0.3) x 10~ 2 , (42) 

bl= 2 = (-0.8 ± 1.2 ± 0.2) x 10~ 2 . (43) 

Our results are compared with previous experiments and theory in Fig. [7J The J x = 2 
results shown are for the two-parameter fit with a 3 = 63 = 0. 



B. Ratios independent of m c and k c 

Using the results from the J x = 2 two-parameter (a 2 ,b 2 ) fit and the J x = 1 fit, we 
find the ratios with the highest statistical sensitivity compare very well with the theoretical 
predictions: 



a 



-£) = 0.67183 = f Cr) = 0-676 ± 0.071 , (44) 

\ Q 2 /CLEO \ tt 2 /th 

'£) = -2.271°- = = -2.27 ± 0.16 , (45) 

,°2 /CLEO \°2 /th 

[ b £] = 0.371^ =( b £) = 1-000 ±0.015, (46) 

V°2 /CLEO V°2 /th 

= -0.111^ = f ^S) = -0.297 ± 0.025 . (47) 

, fl 2 /CLEO V a 2 /th 



C. k c calculation 

Our most sensitive measurement of a magnetic quadrupole amplitude is that of a 2 =1 . 
From theory, we know that 

a 2 J=1 = -^(l±^e) = (l±^)/e, (48) 

where we defined l/£ to be the proportionality between 1 + k c and a 2 =1 . If we use m c = 
(1.5 ± 0.3) GeV, we find f = -(4m c )/E y = -14.0 ± 2.8, so 

1 + k c = £a J 2 =1 = 0.877 ± 0.088 ± 0.034 ± 0.175 , (49) 

where we list the result, the statistical uncertainty, the systematic uncertainty, and the 
theoretical uncertainty from m c = 1.5 ± 0.3 GeV. 
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FIG. 7: Experimental values of the magnetic quadrupole amplitudes from this analysis compared 
with previous experimental values and theoretical expectations. For J x = 2, results are shown for 
the two-parameter fit. CLEO-c results from this analysis are solid circles; Crystal Ball results are 
diamonds [17 L the E760 result is a y 18(, the E835 results are A's 12j, the BESII result is an 
open square 19[], and the theoretical expectations given by Eqs. (f7])- (fTUj) with m c = 1.5 GeV and 
k c = are dashed lines. 
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D. Summary 



We measure significant non-zero magnetic quadrupole amplitudes for the transitions 
Xci —> iJ/^i Xc2 7«Ws an d 4>' — > l'Xci- Our fits to these three amplitudes all agree 
well with the theoretical predictions to first order in the ratio of photon energy to charmed 
quark mass with k c = and m c = 1.5 GeV. The data are consistent with the lattice QCD 
prediction ([TBI) for Xci —> jJ/ip), but not with that for Xc2 lJ I^P [HI- F° r the transition 
if)' — > l'Xc2, we do not measure a significant M2 amplitude, though this case has the largest 
uncertainty since there are fewer J x = 2 signal events and Ey < E-y so I&2I < | ct2 1 - The non- 
zero M2 amplitude in the transitions X(ci,<2) ~^ iJ/tp is evident when comparing the cos# 
histograms for the data with the histograms for phase-space- Monte-Carlo events selected to 
have a pure El distribution and the fitted values of the multipole amplitudes (as shown in 
Fig. [6]). We find that for the J x — 1 and J x = 2 transitions our fitted results differ from the 
pure El value by more than 11a and 6a, respectively. 

The agreement between data and theory is in stark contrast in some cases to previous 
measurements. With about 20 times the largest previous data sample, and a more sophis- 
ticated detector, the matter now seems to be resolved. 
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